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1. INTRODUCTION 

Many optimization problems in practice face obstacles such as the fact that the available data is not 
known exactly when the data problem should be solved. This can be due to measurement errors, modeling 
errors, or a lack of information available when making decisions [1]. Stochastic optimization and robust 
optimization are two methodologies that have been proposed to cope with the issue of data uncertainty in the 
optimization problem. In stochastic optimization, the optimization problem contains expectations and 
probability constraints where the assumptions are a combined probability distribution of known uncertain 
parameters. However, in practice, optimization stochastic does not have sufficient data to estimate its 
distribution, especially in large-scale problems [2]. 

Initially developed in the early 1970 s, robust optimization was intended to give decision makers 
with a framework where probability models could not be worked out [3]. To discover a solution that is viable 
for all potential data parameters up to the worst case, robust optimization assumes that the uncertain data 
have a convex set and a boundary set called the uncertainty set [4], [5]. Specifically, robust optimization 
plays an important role in things that are easy to do and provides an overview of the solution to the worst 
case that efficiently meets the requirements i) high resilience and protection against uncertainty and ii) 
achieving objective values that are close to the problem objectives [6]. One of the well-known methods of 
dealing with data uncertainty in optimization problems is the robust counterpart method proposed by [7]. A 
semi-infinite optimization problem—one with a limited number of variables and an infinite constraint 
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function—is the formulation that is addressed when uncertain data are included in the model under 
consideration [6]. 

Ben-Tal et al. [4], robust formulation can be achieved by assuming the uncertain data set is in three 
sets which are expressed as an uncertain set in the form of boxes, intervals, ellipses, or polyhedral [7]-[9]. 
Ben-Tal and Nemirovski [10], to analyze the computationally robust counterpart analysis, the robust 
counterpart must be converted into linear programming (LP), conic quadratic optimization (CQO), or 
semidefinite optimization (SDO). With this research, the Benders decomposition approach is used to solve 
robust optimization problems for linear programming in a box-interval state in the presence of quasiconvex 
constraints. Where the Benders decomposition method is expected to be able to solve a complex set of 
variables that must be solved in the master problem [11]-[13]. 

— Problem description 

Integer programming problems, such as those posed by Benders decomposition, are well- 
documented in the literature, and it was used to tackle a variety of robust optimization problems, including 
Beale integer programming problems [14]. The Benders decomposition is used by Saito and Murota [15] to 
tackle issues involving robust mix integer programming with ellipsoidal uncertainty. Syed et al. [16] 
considered Benders decomposition to define which of the accessible power plants to reduce price of power. 
Poojari [17], using Benders decomposition to get a robust solution in supply chain planning problems. 
Karamyar [18], using Benders decomposition for location problems and machine allocation in health care 
schedules. Emami [19], developed a robust counterpart model in the problem of receiving orders and 
scheduling in machines. Order issues are determined by different date, revenue, late penalty, processing time 
and setup time on the machine. From the proposed model, several valid pieces are introduced to accelerate 
the convergence of the Benders algorithm and heuristic methods to get a feasible solution. 

From the various literatures described above, this study discusses solving optimization problems for 
quasiconvex constraints in a box-interval uncertainty. In this problem, robust counterpart is expressed as a 
semi-infinite linear programming problem which is an optimization problem with many variables and 
infinitely many constraints [20]. This method will create a strong companion model of the proposed model 
by introducing several valid parts to accelerate the convergence of the classical Benders algorithm and 
heuristic methods to get a feasible solution. 


2. THE PROPOSED METHOD 
2.1. Robust optimization 
The formula for the indefinite linear optimization is as follows: 


min{c"x: Ax < d}ccaayeu o) 


where c € R”, A € R™*" and d € R” is an uncertain coefficient, and U shows the user the uncertainty set. 
The basic paradigm of robust optimization is based on the following three assumptions [3], [4], [21]: i) all 
decision variables x € R” represents the here and now decision, where the decision variable must get a 
certain numerical value as a result of solving the problem before the actual data, ii) it is the responsibility of 
the decision maker, if the real data falls within a predetermined uncertainty range U, to make a decision. Iii) 
the constraint of the uncertainty problem is hard, i.e., the decision maker cannot abide the violation of the 
constraint when the data is in a predetermined uncertainty set U. 


2.2. Robust counterpart 
Assume c E€ R” and d E R” is certified, then the robust reformulation of (1) is called the robust 


counterpart problem which as shown in [4], [22]: 


min{c?x: C(x < q,Y E Z} (2) 


where Z c R! shows the primitive uncertainty set. A solution x € R” is called robust feasible if it satisfies 
the uncertain constraint C(¢)x < q for all realization ¢ E Z. 


2.3. Box-interval uncertainty 
In box-interval uncertainty U, an optimization problem can be formulated as a semi-infinite 


program with robust counterparts [23]: 


minz (3) 
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s.t.f (x,0) < z,YÇ E€ U, 
g(x,f) < 0,YÇ € U, 
xEX 


In (3), there is Us» = {[]|@|].. < Y} [6] where w is a control parameter in the uncertainty set [24] as 
illustrated in Figure 1. Based on [25] the value of Y < 1. 


Classical Set Flexible Set Cut Set 


Figure 1. Box-interval illustration 


2.4. Convergence in robust optimization 
Definition 1: convex function: a function g(p, ñ) is called quasiconvex on p E [—Ap,+Ap] if and 


only if 6 E [—Ap, + Ap], g(p, ô) < max{g(p, Ap), g(p, —Ap)} for all p [20], [26]. 
Definition 2: objective robustness: n the case of a candidate points (p°,q°) objective robustness 
predominates if the inequality is lowered: 


CCRA Cc Ac 
f(p%a pore ,q°,0,0) <1 (4) 


fulfill for all 6 € [—Ap, Ap] and ĝ E [—Aq, Aq]. 
Definition 3: Robustness of feasibility: In the case of a candidate solution (p‘°, q°) if: 


g; qp, G) < 0, Vj = 1, ...,J (5) 


fulfill for all 6 E€ [—Ap, Ap] and ĝ E [—Aq, Aq], then the robust feasibility is prevailed. 
Assume that the objective function is convex with the constraint convex so that the solution persists 
when nominal p = g = 0. 


min f (p,q) (6) 
pq 


s.t. 00.0019) < Oj = 1, J 
p E R”,q E Z™, p E€ R", ĝ E€ R™ 
vô € [-Ap, Aq], VG € [—Aq, Aq] 


Definition 4: robust points: points that meet robust eligibility for (6) are considered robust points. Sp 
represents the set of robust points for (6). 

Definition 5: locally optimal robust: for the robust optimization problem, the locally optimal robust 
p * solution is a robust point such that there is a set of neighbors U of the optimal p * robust solution 
f(p *) Sf (p), Vp € U. 

Definition 6: globally optimal robust: for formulation (6), globally optimally robust (p*,q*) is an 
optimal robust point on f (p*, q*) < f (p,q), V@, gq) E Sr. 


3. METHOD 
3.1. Robust optimization re-formulation 

By adding an auxiliary function a, (6) will produce two master problems for the defined Benders 
decomposition: 
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ay (p,q) = mnp (p,q) (7) 


8.t.9,0,¢3,9) < 0,j = 1, J 
p E€ R”,q E Z™ 


a_(q) = minf (p,q) (8) 
p 


s.t. gjp,q ĝ,ĝî) < 0,j = 1, J 

pe Prixea ĝ = fixed 

p ER” 
3.2. Algorithm improvement 

By following the basic algorithm of [27]-[29], the following Benders decomposition improvement 
algorithm is used to solve problems (7) and (8): 
Step 0. Set the / iteration to 0 with the smallest constant as the tolerance value. 


Step 1. Set iteration for I = ite +1 where variables ($,q) are a complex variable that is fixed in the 
subproblem. So, the master problem becomes: 


mina, + 
ane faa) 


s.t. Adown < q < dup 


aq < amin (9) 


—Ap <q < Aq 
aq Say (10) 
the limits of @p, @q determined and given a, = a#°, aq = ag? and q = qifen P = Pikea 4 = Yared: 
Step 2. Fix complicated variable values p, then solve the subproblem with Bender's standard decomposition 
method. 
w= minfp(p) (1) 


s.t.9,0,¢3,9) < 0,j = 1, J 


4 = ixea (dual = 4°) 


A — fite 
p= Prixed 
a _ aite 
q = a fixed 


Step 3. Check for convergence. Set Zub =W, Zmasti = Oy Zales = a If the difference z 
Zdown1 Š € and Zup — Zdownz < € then the iteration stops. 
Step 4. Add standard Benders snippet for master problem (9). 


up ` 


a 2h (pie ae BRL GER) + a’ (a — af (12) 
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Step 1 (return). N represents the number of iterations that have been completed. Once Benders pieces have 
been added, solve the following master problem, 


minaq T fal). t Adown < q = dup (13) 
q” 


l l asol asol it l 
ap = f (PR, qil, PRL, aie’) + 2° (q — qi 
for ite = 1,...,N 


min 
Qq Z Qq 


q 


return to step 2 and continue with iterations until convergence is met. 


4. RESULTS AND DISCUSSION 
Numerical calculations are show the improvement algorithm for the Benders decomposition method 


to obtain a robust solution to problems (7) and (8). In this case, the robust optimization problem lies in the 
box-interval uncertain constraint. 


minz = —=p-—q (14) 
Pq 4 


s.t.(-1+%,)p+(-1+@%)q <5 
(-2+%)p ++ Gq << 


(F+6)pt+U+isqs= 

(1 +p + (—1+¢)q < 10 

(1 +)p < 16 

(1+ %)q 20 

vč; E box — interval,i = 1, ...,I 


4.1. Solution 


By taking ¢ € [—0.1,+0.1] and ¢; = 0.1 which is substituted into the problem constraint (14), 
problem (14) will become a deterministic problem (15) which is robust by forming a robust counterpart into 
the set of box-interval uncertainties: 


minz=—7p—4q (15) 
s.t. —0.9p — 0.9q <5 

—0.4p + 1.1q < 7.5 

0.6p + 1.1q < 17.5 

1.1p — 0.9q < 10 

1.1p < 16 


1.1q 20 


the next stage, the deterministic optimization problem (15) is solved by the improvement algorithm Benders 
Decomposition starting with: 
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Step 0: iteration v = 1. The master problem, 
TRS = =p +a 
s.t.1.1p < 16 
—25<a 


Step 1: the subproblem solution: 


min z = —q 
q 


s.t. — 0.9p — 0.9q < 5 

—0.4p + 1.1q < 7.5 

0.6p + 1.1q < 17.5 

1.1p — 0.9q < 10 

1.1p < 16 (dual = A) 

11q>0 

The solution is SY = (p®, q) = (14.5,8) by minimizing the objective function z. And the 


optimal value of the dual variable with the constraint pY = 14 is AY = 0.5. 
Step 2: check for convergence. The upper and lower limits of the optimal objective function, 


Zap = — {pO — q@ = —2 (14.5) — 8 = -11.625 


Zom = — PD +a) = —2 (14.5) — 25 = -28.625 


because of the difference cae E zD n = 17 > g, then the process continues to the next step. 


Step 3: master problem solution. Update the iteration v = 1+ 1 = 2. Master problem: 
minimize = <P +a 
s.t.8.7 + 1.1(p — 14) < 17.5 
1.1p, < 16 
& < 17.5 


where the first constraint is done by Benders cut 1 which is related to the previous iteration. 

Because too many iterations have been carried out and the convergence has not been achieved, the 
problem solving (15) is continued with the help of the production and operations management—quantitative 
methods (POM-QM) software. By using the POM-QM software, the convergence results are obtained with 
the upper and lower limits of the optimal objective function being: 


2) = -1p® — q® = -12.95 
ONPE -1p® +a® = —12.95 


Zdown 
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because of the difference z® = zE n = —12.95 + 12.95 = 0 < € then the optimal solution is obtained by 
p *= 10,q *= 10.45,z x= —12.95. The results of the iteration summary with the improvement of the 
Benders decomposition method are shown in Table 1. Figures 2(a) and 2(b) provide examples of the 
problem's feasible and robust feasible regions problem (14). 


Table 1. Optimal solution and iteration improvement of Benders decomposition 


Iteration p® q® a® 24 z 2 a 
1 14.5 8 -25 0.5 -11.625 -28.625 
2 0 4.5 -17.5 -1 -4.5 -17.5 
3 4.5 8.5 -13.33 -0.5 -9.625 -14.455 
4 -6 0 -5.25 -0.5 -1.5 -6.75 
5 10 10.45 -10.45 -0.5 -12.95 -12.95 


(a) 


(b) 


Figure 2. Provide examples of the (a) illustration of feasible region and (b) illustration of robust feasible 
region 


The deterministic problem (14) is a linear programming problem by giving the value of 
p *= 10,q *= 12.5 and z = —15. It is possible to establish the robust optimization (14) solution by 
examining the corner point of the robust workable area, which provides the value of 
p *= 10, q *= 10.45 and z += —12.95. Thus, as illustrated in Table 2, the globally optimal interval is met. 
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Table 2. Comparison of feasible regions with robust feasible regions 
Description Deterministic solution Robust solution 
p* 10 10 
q* 12.5 10.45 
zZ* -15 -12.95 
CONCLUSION 


To create a stochastic optimization model, one assumes that the input data is known precisely and is 


equivalent to a predetermined nominal value. Data uncertainty has a significant influence on the model's 
quality and feasibility, which this technique does not consider. When data are obtained from various nominal 
values, certain constraints may be broken and the best solution established using nominal data may no longer 
be optimal or feasible. Problems that include parameters or decision variables with box-interval uncertainty 
may be solved using this method. It is suggested that Bender's decomposition improvement approach be used 
to produce a robust optimal linear programming solution. 
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